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We explore in detail the semiclassical environment of collapsing shells of matter, and determine the 
semiclassical flux measured by a variety of observers. This study is a preliminary step in a broader 
investigation of thermodynamic properties of the geometry of collapsing objects. Specifically, in 
this paper we consider spherically symmetric null and timelike collapsing shells which form an event 
horizon, and calculate the flux measured by observers both inside and outside the shell, and both 
inside and outside the event horizon, and find nontrivial results in most of the cases. Additionally, 
^ ' we also investigate the environment of a shell which collapses but does not form a horizon, halting 

O ' at some radius larger than the Schwarzschild radius, and find that such an object generically gives 

, rise to a pulse of radiation which is sharply peaked as it travels inwards and is reflected at the 

origin, and eventually emerges from the shell in a "thermalized" form. Our results have potential 
(N ■ consequences in addressing questions pertaining, e.g. to black hole entropy and backreaction. 
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I. INTRODUCTION 



The behaviour of quantum fields in the spacetime geometry of black holes and gravitationally collapsing objects in 
general has been a subject of detailed investigation since Hawking's seminal discovery [l| of a nearly thermal flux 
of particles that is measured by observers located outside the collapsing object. Several authors have subsequently 
contributed towards building an understanding of the semiclassical environment of collapsing objects 0, (sj, y, 
\ One of the most fascinating aspects of this field of research is the fact that classical black holes appear to behave 
■ like macroscopic thermodynamical objects An important offshoot of this behaviour is the so called information 
[ paradox (see, e.g. for a recent interesting discussion), and an understanding of this thermodynamics at the quantum 
• level is also likely to have consequences towards building a quantum theory of gravity 8]. 

There are several questions that can be asked within the semiclassical framework. For example, considering a 
0^ , Schwarzschild black hole for simplicity, it is a celebrated fact that the black hole can be attributed an entropy which 
' scales as one quarter of the area of the event horizon in Planck units. This entropy however has a physical relevance 
J> ] only to those observers who choose to remain outside r = 2M forever, i.e. - for observers who truly experience r = 2M 
as an event horizon. Since forever is a long time to wait, it becomes interesting to ask whether the concepts of black 
hole entropy etc. can be meaningfully addressed locally. For example, what entropy (if any) would be attributed to 
I the surface r = 2M by observers who remain outside for an arbitrarily long time, before suddenly choosing to jump 
■ " ' into the black hole? Such observers will clearly eventually have access to many more microstates compared to their 
cousins who choose to remain outside forever, but until the decision to fall in is actually made, there are no grounds 
to distinguish between these two classes. 

Returning to the context of collapsing objects, further complications arise when one begins to account for the 
backreaction of the semiclassical fltrx of particles, on the classical geometry of the collapsing object. The basic idea 
here is that the loss of energy from the collapsing object (which we take to be spherically symmetric, for simplicity), 
in the form of the semiclassical flux, implies that the correct background geometry in the exterior of the object must 
be time dependent, rather than the static Schwarzschild spacetime 0, H, • Due to technical challenges, the exact 
effect of the backreaction has been a subject of much debate, with claims to the extent that the formation of the 
event horizon itself may be hindered by the backreaction [ll|, • While this appears unlikely [l^ , it is nevertheless 
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true that determining the exact geometry in the late stages of black hole evaporation, remains an open problem. 
Consequently, questions regarding the thermodynamic behaviour of the geometry of the black hole/collapsing object 
also remain clouded by the backreaction. 

In this paper, as a preliminary step towards addressing the broader issues raised above, we set a simpler mandate 
for ourselves. We wish to explore the full semiclassical environment of a collapsing object, which we will consider to 
be a spherically symmetric thin shell [2| • Classically, the exterior of the shell is described by Schwarzschild spacetime 
while the interior is Minkowski spacetime. In particular we wish to determine the semiclassical flux, as encoded in 
the renormalized stress tensor, observed by a variety of observers, both inside and outside the shell, before and after 
horizon formation. 

There are several motivating factors behind such a study : Firstly, while the flux seen by observers outside the 
shell and the event horizon has been studied in great detail (beginning with Hawking's own calculation), the nature 
of the flux inside the shell, and inside the event horizon but outside the shell, remains largely unexplored. In fact, an 
argument one usually encounters is that observers who fall into a black hole do not encounter a significant flux. Our 
calculation assumes zero backreaction, and we will see that in this case an infalling observer will see a flux of outgoing 
particles which in fact increases without bound as the observer approaches the singularity (although as expected it 
is finite at all other locations, including r = 2M). Further, at least in the case of timelike collapse, we find that 
observers inside the shell also see a flux of particles, which might have some interesting consequences which we will 
return to discuss in the final section. 

Apart from collapsing shells which form an event horizon, we will also discuss the case, first described in Ref. 13?|, 
of a shell which even classically does not form a horizon; instead it asymptotically approaches a final radius larger 
than 2M. In Ref. such a shell was shown to have very interesting properties, although that discussion was limited 
partly to order of magnitude estimates, since some required expressions were not analytically obtainable. In this 
paper we will perform a full numerical study of such a shell trajectory, and also explore its semiclassical behaviour 
both inside and outside the shell. We will find that such trajectories are generically expected to generate a pulse 
of radiation which travels inwards, reflects at the center and travels back outwards, eventually reaching an outside 
observer in a "thermalized" form (to be made precise later). 

At this stage, it is worth addressing why we choose to study the seemingly unrealistic situation of a collapsing shell. 
The first reason of course is that the geometry of a shell is simple to handle analytically, not least because the interior 
is just Minkowski spacetime. A second important reason is that the absence of interior matter allows us to study 
the purely gravitational aspects of quantum fields in the interior. It is important to understand such effects before 
addressing the problem of quantum fields whose semiclassical effects interact with those of classical matter inside a 
star for example. Throughout our discussion we will simply prescribe the classical trajectory of the shell, and not 
worry about the specific (singular) stress-tensor needed to generate the trajectory. This is not expected to interfere 
with the conclusions we draw regarding the semiclassical behaviour of the shell. 

The plan of the paper is as follows : In Section[TT]we start with the example of a null shell which has been collapsing 
forever. While such a geometry has been studied in the literature we note that the region of Schwarzschild 

spacetime for r < 2M has been largely ignored. Another reason to study the null shell is that all expressions for the 
flux, for all observers, can be obtained analytically, with no numerical inversions involved. This section will therefore 
serve as a useful warmup to the more involved calculations for timelike trajectories. In Section [TTTl we turn to studying 
timelike trajectories, where many of the calculations will be numerical. We will first study a trajectory which forms 
a horizon, describing the flux seen by various observers, and then turn to the trajectory mentioned earlier, which 
asymptotically halts without forming a horizon. We will end in Section IIVI with a brief discussion and prospects for 
future work. 



A. Notation 

To avoid confusion, in this subsection we will lay down all the notation that will be used subsequently and in the 
next subsection we recapitulate the prescription for calculating the semiclassical flux in a given geometry. Note that 
all our calculations assume zero backreaction. 

We will be assuming Planck units throughout, with G, h and c set to unity. We will mostly follow the notation 
used in Ref. p^ . which in turn was based on the notation in the review by Brout et al. [9i|. Additionally, to avoid 
cumbersome factors of 2M (where M is the mass of the collapsing object), we will assume that all the coordinates 
are non-dimensionalized by dividing out by 2M. This will lead to expressions for the flux, etc. with factors of 2M 
missing, which can be put back appropriately by dimensional arguments. Latin tensor indices a, 6, .. range over 0..3, 
and we use the mostly plus sign convention. 

As mentioned earlier, the exterior of the shell is Schwarzschild spacetime while the interior is Minkowski, and we 
can use the same radial coordinate r to foliate both regions. We refer to Schwarzschild time as t and Minkowski time 
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as T. Outgoing null coordinates will be a variation of the letter u, similarly ingoing coordinates and v. To be precise, 
the interior Minkowski null coordinates will be 

V = T + r ; U = T-r, (1) 

and the exterior Eddington null coordinates will be 

v — t + Xf,{r) ; w = t — x,(r), (2) 

where we have defined the tortoise function 

x,(r) =r + ln|r- 1|, (3) 

(which should be read as a:*(r)/2M = r/2M + ln |r/2M — 1|). All "Kruskalized" (i.e. - exponentiated) versions of the 
null coordinates will be referred to using calligraphic fonts. Specifically, lA and V will denote the standard Kruskal 
null coordinates 

V = 2e''/2 ; W^TSe-''/^ ; ^YV 4p7(1 - r) , (4) 

with the upper sign in U for r > 1 (the Schwarzschild radius being unity in non-dimensionalized form) and the lower 
sign for r < 1. Other Kruskalized coordinates will be prescribed as needed. 

With this setup, the (non-dimensionalized) metric in various regions and in various coordinates becomes (suppressing 
the line element for the 2-sphere, which is always r^dil^), 

dsf^j ^ -dT^ + dr^ = -dUdV , 



ds^^t = -(1 - l/r)dt^ + (1 - l/r)-^dr^ 

= - ^1 - dudv = --^dUdV . (5) 

The subscript s will always refer to evaluation on the shell, so the trajectory is r — Rs, v — Vs, etc. The specific 
parameter used along the trajectory will be chosen according to convenience, and could be one of the coordinates. 



B. Recap of semiclassical prescription 

Throughout, we will be interested in the semiclassical flux J- — na{T^)u^ measured by an observer with 4- velocity 
M° where Uau"' = 0, Uau"' = 1 = —Uau", and (Tab) is the renormalized stress-tensor vacuum expectation value 
(VEV) of a quantized, massless scalar field in a specific vacuum state. We refer the reader to Refs. Ky, ^ for detailed 
treatments of the formalism, and to Ref. [isj for a pedagogical introduction. The Hawking flux has the value 
(7r/12)r^ = (1927r)^^(2M)~^, where Th — I/SttM is the Hawking temperature. The vacuum state |in) to be used, 
is completely defined by specifying its behaviour at past null infinity X~ and by requiring that its mode functions 
vanish at reflection at the origin r = 0. Typically one requires that the vacuum on I~ be identical to the Minkowski 
vacuum, so that the ingoing modes of the state |in) become ~ e"*'^'" for a typical timelike collapse. We will see that 
the situation for a null shell collapsing forever, is somewhat different, since in this case the portion of T~ relevant for 
thermal flux in the exterior, lies inside the shell. However, here also we will be able to meaningfully compute the flirx 
observed by specific classes of observers. 

An important assumption we make is to work in the 1-1-1 s-wave approximation, ignoring the Schwarzschild 
potential barrier, so that the results of 2-dimensional conformal field theory (CFT) apply. For a discussion of this 
approximation and its validity, see the review by Brout et al. Q. In this approximation, for a vacuum state |vac} 
with ingoing and outgoing modes given respectively by e~^^^^ and e~'^^^+ , such that the metric in the relevant region 
in terms of these null functions is ds^ = —C{f+,f-)df-df+, the components of the renormalized stress-tensor VEV 
in the state |vac) are calculated as 

(T/_/_)vac = -F/_[C(/_,/+)] ; {Tf^uU,^-Ff^[C{f^,f+)], (6) 

where the derivative operator F is defined, for any function f{x), as 

F,[f{x)]^ ^f'/^dlr"^ ^ -^bdllnf - {d.Anff] . (7) 

izTT 4o7r L J 
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FIG. 1: Penrose diagram for the null shell that has been collapsing for all time Ingoing rays have V = const and outgoing 
rays have U = const. 



The cross terms {Tf_f^ )vac can be determined by requiring that the stress tensor be covariantly conserved. The 
following relations are useful while analyzing the flux in different regions : 

FAx]^-^, (8a) 

F,[l/fix)] = -F^lfix)] + ^{d^lnff , (8b) 

FMix)9i^)] = F,[/(x)] +F,[5(x)] + ^(9, In/) (9, In 5) , (8c) 

Fs[fix{s))] = ^ (f,[/(x)] + -1- (a. In/) (a. Ins')) , (8d) 

where in the last equation, s' = ds/dx. 

Hawking's result follows from recognizing that for a generic timelike collapse, requiring the vacuum |in) to be 
Minkowskian at 1~ , the corresponding modes in the exterior region become e"*"^", e~'"'^("\ where G{u) is that value 
V = G for which an ingoing ray of constant G, becomes a ray of constant u after being reflected at the origin and 
emerging from the shell. Further, one finds that due to the presence of the event horizon, the generic behaviour of 
G{u) for M ^ cx) is dG/du ~ e~"/^, which leads to the result that in the state |in), an exterior observer using the 
Eddington u,v coordinates measures a constant outgoing thermal flux (T„„ ) — {Tyy ) = (7r/12)T^ as u ^ co. [The 
notation ( Tab ) with no subscript will always refer to the vacuum state |in).] We will see this result in the trajectories 
we study. For further details see Refs. [sl. M. [Tsj . 



II. NULL SHELL 

We begin by considering a null shell that has been collapsing for all time [1, [3, [iBl ■ The global structure of this system 
is different from the usual timelike collapse situations which we will discuss in Section IIIIl since a part of in this 
case lies inside the shell. Since the matching of metrics in ([5]) occurs along a line of constant V (which we conveniently 
choose to be y = 0), there is a freedom in relating the exterior and interior ingoing coordinates. We exploit this 
freedom to set Vm = 2e^/^ ~ V. Additionally, we use the Kruskalized outgoing coordinate Z//„i = 2e~'^/^(l + U/2), 
and it is then easy to show that on the shell (i.e. - at Vm = V = 2) we have Um = U for all r. In other words, we can 
use the usual Kruskal coordinates lA, V over the entire spacetime, with the metric given by 

dsL = -dUdV = ~^^^dUdV , 

dsl^,^^'-^dUdV, (9) 

with the event horizon r = 1 corresponding to the outgoing ray lA = Q ov U = —2. Figure [T] shows the Penrose 
diagram for this spacetime [14| . 
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U/2M 



FIG. 2: The semiclassical flux (normalized by the Hawking value) measured by static observers outside the event horizon, for 
the null shell trajectory. 



In the usual timelike collapse situations where the shell would typically start collapsing from some finite radius, X~ 
lies entirely outside the shell and the spacetime is asymptotically flat. It is then natural to choose a vacuum whose 
ingoing modes at are ~ e^*"", v being the ingoing Eddington coordinate. In the present situation, it is natural 
to define the vacuum in the interior region of X~ using the Minkowski ingoing coordinate V. Note however that our 
choice of mapping the ingoing coordinates between the interior and exterior, ensures that V = v at the boundary, 
and hence the |in) vacuum can still be chosen to have modes ^ g"*'^'" — e~*'^^ on all of 2^ . The reflection condition 
at r = then fixes the outgoing modes of |in) to be ^ g-'-'^u ^j^g interior region. Equivalently, the function G 
discussed in Section ITbI is simply G{U) = U in the interior. It then follows that the renormalized stress tensor in the 
state I in) is identically zero in the interior of the shell, since the conformal factor in the interior is constant in the 
U,V coordinates. 

In the exterior, for r > 1 we want the flux measured by static observers Oi who use the Eddington u,v coordinates. 
The function G{u) is now determined by the implicit equation 



The stress tensor components in the state |in) are given by 



-^2 \ 1 / . 3 
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dG_ 

du 



(10) 



(11) 



[Note that the r-dependent part of these expressions is simply the static contribution in the so called Boulware state.] 
The outgoing flux measured by observers Oi is 



dG 



du 



1 



1 



12ttG{u) 



(3 + 2G(w)) . 



(12) 



We can analytically see that as m ^ oo, G ^ —2, dG/du ~ e~"/^, and hence the flux asymptotically approaches 
the Hawking value (note that in our convention for non-dimensionalization, the Hawking value is l/1927r). Figure [2] 
shows this flux (normalized by the Hawking value) as a function of u. 

An interesting region to consider is the exterior of the shell for r < 1. In this region it is more natural to use the 
Kruskal coordinates U,V, since this region is likely to be accessed by freely falling observers O2. The stress tensor (in 
the entire exterior region r > Rs) in these coordinates is 
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[Tvv ) - 7^(r™ ) , 



[ Tuu ) = Tj^i Tuu ) 



(13) 



Although superficially divergent, the expression for {Tuu ) is actually finite on the horizon U = 0. The structure of 
( Tyy ) and ( T„„ ) shows, on the other hand, that observers who fall into the singularity will generically see a flux 
that increases without bound. We will see this feature again when studying timelike trajectories in Section [TTTl below. 
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where we will display a more detailed expression and a plot for the flux seen by observers 02- A third class O3 which 
remains in free fall inside the shell for a certain time, does not see a flux in the present null case. The timelike case 
however will prove to be different. 

We end this section by noting that although inertial observers in the interior of the null shell do not see a semiclassical 
flux, there is a natural class of accelerated observers who do see a flux. These are observers following trajectories 
picked out by the Kruskal null coordinates, and are related to Rindler observers in the Minkowski region. (Notice 
the contrast with the exterior, where such Kruskal observers are inertial, while observers following trajectories picked 
out by Eddington coordinates are Rindler- like.) For a vacuum state |vac ) defined by the modes e~*"^, e"''^'^ in the 
interior, the stress-tensor components are given by 

(7vv)vac = (^2"^^) ' 

{TuuUc^[^Tjj)^{U^-W+l2) . (14) 
III. TIMELIKE SHELL 

We now turn to the more realistic situation of a shell collapsing along a timelike trajectory. We will discuss two types 
of trajectories : the first is a "conventional" trajectory in which the shell forms a horizon in finite proper time, and 
continues to collapse to a singularity. The second is an "asymptotic" trajectory, in which the shell asymptotically 
approaches a final radius larger than its Schwarzschild radius, so that a horizon is never formed. In Ref. [l^ the latter 
was shown to have very interesting behaviour, in that an outside observer would see a constant Hawking-like flux for 
some finite time interval, after which the flux would exponentially decay to zero. In Ref. ^13 ] however, only a part 
of the flux was explicitly calculated, with order of magnitude estimates for the remaining part, since the functional 
inversions needed could not be performed analytically. Below we will show the results of a full numerical calculation 
of the flux, which will corroborate the results of Ref. JL3] ■ We will also see some additional effects as mentioned in 
the introduction, which are only revealed in the numerical analysis. All numerical calculations were performed on 
Mathematica. We begin with a discussion of the conventional horizon-forming trajectory. 

A. Horizon-forming trajectory 

We parametrize the trajectory using the shell proper time. Our trajectory is r = Rs{t) (assuming a non- 
dimcnsionalised form for all coordinates), where 

R's=m^ ; Rs{r) = Ro, r<0, (15) 

where a prime denotes a derivative with respect to the argument t, 9 is the Heaviside step function, and / is a function 
which is chosen so that / ^ I at late times. Note that R'^ < 0. The trajectory therefore remains at a constant radius 
Rq until T — 0, and eventually becomes a pseudo-free-fall curve, mimicking a particle freely falling in the potential of 
a body of mass M. The specific form of / we choose is 

fir) = 1 - e~(^^)" , (16) 

with K and a as parameters. We will fix a = 5, which ensures that Rs{t) has a continuous third derivative at r = 
(and hence for all r < Tmax, where i?s(Tmax) — 0). Since the parametrization is in terms of the shell's proper time, the 
equations (O are sufficient to completely determine the trajectory in terms of other coordinates u, v, etc. as well. In 
particular, the functions Vg, Ug are useful in the analysis, and can be obtained by integrating the following expressions 

v', = llg ■ g{T) = -Ri + ^Rf + 1 - l/R, , 

c/: = -i?:, + /RfTT, (17) 

and when needed, we have Vg — Ug -\- 2Rs, Ug = Vg — 2x^{Rg). 

As in the null case, we will calculate the semiclassical flux in this spacetime, for three classes of observers Oi, O2 
and O3 . Observers Oi and O2 always remain outside the shell, while observers ©3 remain inside the shell until some 
time. More precisely, Oi correspond to r = const observers with r > Rq, O2 are free-fall observers who cross the event 
horizon and fall into the singularity, and O3 are chosen to be r = const observers inside the shell, who are tracked 
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FIG. 3: Penrose diagram for timelike collapse of a shell, with Ro = 6M, 2Mk, — 0.3. The origin (leftmost thick black curve) 
remains regular until the shell collapses to the singularity (thick dashed black) . Various lines correspond to the shell trajectory 
(thick solid red), the surface r = 2M (thin solid blue), the last escaper null ray (long dashed green), the ray which enters at 
r = (short dashed green), and example trajectories of observers Oi (r = 6.15M, exterior short dashed black), O2 (/? = —0.5, 
7 = AQAM, fine dashed purple) and O3 (r = M, interior short dashed black). 



until the time when the shell crosses their location. Figure [3] shows the Penrose diagram for this spacetime, with 
example trajectories of Oi, O2 and O3 observers also shown. Recall from Section ITB] that the modes corresponding to 
the vacuum state |in) for this spacetime, are completely fixed by requiring that (a) the modes on X~ become ~ e~*'^", 
and (b) the modes vanish at the origin r = 0. Given these conditions, we can calculate the function G, which then 
determines the form of the conformal factor used to calculate the flux in Eqn. 



Flux for observers Oi 



The calculation in this subsection is well known in the literature (see, e.g. Ref. [31). We will nevertheless display 
some relevant details which demonstrate the general procedure followed in the subsequent calculations. For observers 
Oi with constant r who use the Eddington u as their natural retarded time, the function G(u) can be parametrized 
by r, with u = Us(t), and is given by 



(18) 



i.e., a ray which exits the shell at r, must have entered at a time ti determined by the second equation in (|18p . which 
is just the reflection condition. G is then by definition the value of v at the entry point ri, and is treated as a function 
of T. Further note that 



dG 
du 



dvsiri) dri 



1 



dri dr (dus/dr) 



= T^'s(n)/<, 



(19) 



where we have introduced the convenient shorthand notation of using a prime to denote a derivative with respect to 
the argument, and evaluation of all functions will be at their natural argument r, unless explicitly stated as in the 
case of v'g (ti ) . We will use this type of analysis for the observers O2 and as well. 

The (outgoing) flux T for these observers is independent of r, and is given by (see also Ref. Q) 



Tiu) ^ -{TI) ^ {Tun) - {T,v) ^ F.^ 



dG 



du 



(20) 
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U/2M 



FIG. 4: The flux (normalized by the Hawking value (7r/12)r^) as a function of it, for the timelike horizon- forming trajectory 
with Ro — 8M, for r = const observers outside the shell. The two curves correspond to 2Mn — 0.1 (solid blue) and 2Mk — 0.3 
(dashed purple). The initial behaviour of the flux depends on the trajectory details, but eventually the flux attains the constant 
Hawking value (short dashed black line). 

where the derivative operator on the extreme right was defined in Eqn. ([7]). Using the relations in Eqns. ([5]), and 
parametrizing again using r, we can rewrite this as 

Hr) = ^\^-FAu',{r)]+T[^Fy[v'M])\y=r,ir)+Fr[r[ir)]'^ . (21) 

The expression on the right can be evaluated numerically, with the only nontrivial part being the inversion involved 
in determining ti(t). The result is shown in Fig. 31 where the flux is plotted as a function of u. As expected, this 
flux attains a constant value (7r/12)T^ at late times. 



2. Flux for observers O2 

We take observers O2 to be freely falling into the horizon and eventually into the singularity. Defining the timelike 
and spacelike Kruskal coordinates 

tK = l{V + U) ; rK = ^(V-W), (22) 
such an observer can be described as moving at constant velocity in the (tx, ''k) coordinates, so that 

rK = PtK + 7/2 , (23) 

for some constant 7 and /3 < 1. This trajectory can also be parametrized using the shell's proper time r, and we 
ensure that this observer always remains outside the shell by choosing an appropriate 7 and tracking the trajectory 
only after a fixed time Tob.min- (Recall that for the eternal Schwarzschild black hole, a trajectory such as Eqn. ((23l) 
corresponds to an observer who pops out of a white hole at t = —00, reaches a maximum radius r = r„iax, and 
falls into the horizon at t = +00.) It is convenient to write the trajectory in terms of the r and V coordinates, as 
r = -Rob(T), V = Vob(T), where 

Vob(T) = ((l+m+7)/(l-/3), 
4e«°Hl - i?ob) = i^sVob , (24) 

where we have used the relation between r and the Kruskal Uy given by lAV = 4e'"(l — r), and Us{t) is the trajectory 
equation which can be obtained from the matching conditions on the shell, 

Us{T)=2e''^-^^'\l-Rs). (25) 

The outgoing flux F for such an observer becomes 

F^ na{Tl)u\^^ = - (re'-{(rvv) - {Tuu)])U . (26) 



Flux 



9 




FIG. 5: The normalized flux (red solid line) for the timelike horizon-forming trajectory with 7?o = lOM and 2Mk — 0.1, for 
a freely falling observer with /3 = —0.75 and j = 27. 7M. The flux is plotted as a function of retarded Kruskal time U, with 
horizon formation occurring at W = 0. The vertical dashed line indicates the time at which this observer hits the singularity 



SO that the flux normahzed by the Hawking value, denoted J^, is 

^^-[^(^^)-^-o.^^)]L--(-«[i)L. 

Note that while it seems analytically easier to parametrize this trajectory using the Kruskal coordinate U, numerically 
it is more convenient to use the shell proper time, due to the appearance of the function G in the flux. Since we are 
parametrizing the observer's trajectory using the shell's proper time, r and V in Eqn. ([27)) will be replaced by Roh{T) 
and Vob(T) respectively, and the last term involving the F operator has a form identical to the RHS of Eqn. (j21|) . 
with Us(t) replaced with Us{t). 

This flux is also straightforward to compute numerically, and the results are shown in Fig. [51 for a shell trajectory 
with i?o = lOAf and 2Mk = 0.1, for an observer with (3 = —0.75 and 7 = 27. 7M. We see that while the flux is small 
(with magnitude comparable to the Hawking flux) as long as the observer remains outside the horizon, and remains 
finite at horizon crossing as well, once the observer approaches the singularity, the flux increases without bound. This 
is also obvious from the expression (|27p . where the first term in square brackets is singular as r ^ 0. 



3. Flux for observers O3, 



We now turn to the interior Minkowski region, and consider the fiux measured by r = const observers during the 
time that they remain inside the shell. We do not track such observers once the shell has crossed them, since this 
event is expected to be singular for such observers even classically, and we have not been accounting for the classical 
stress-tensor of the shell in our analysis. Unlike the null case, here we will see that interior observers do detect a 
nonzero fiux, which is a consequence of the difference in global structure of I~ in the two cases. 

The r = const observers inside the shell use r = (y-|-[/)/2 as their natural proper time. For any event T' = (C/*,14) 
along the trajectory of such an observer, there are two events Vi and V2 relevant for the flux calculation. Vi is defined 
as the entry point (into the shell) of a ray which, after reflection at the origin, becomes the outgoing ray U — U^. F2 
is similarly the entry point of the ingoing ray which after entering the shell becomes V = V^,. Let ti and T2 be the 
values of the shell proper time corresponding to the events Vi and respectively. In terms of (T, r) we have 

V,iri)^U,^F~r ; V,{t2) ^ V, = T + r . (28) 

The modes of the |in) vacuum are then given by ~ g-i^vA'r2) _g-jwtis(Ti)_ Using the CFT rules for calculating the flux, 
it is then not hard to show that the normalized outgoing flux T = {{Fiju ) ^ {Tvv )) / i'^T'^ / 12) for such observers is 
given by 

.F = 5(ti)-5(t2), (29) 

where the function S{t) is defined as 

1927r 

= (FrHir)] - F.K(r)]) . (30) 
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FIG. 6: The normalized flux for the timehke horizon-forming trajectory with Rq — lOAf and 2Mk — 0.1, for three r = const 
observers inside the shell, as a function of Minkowski time T. The curves correspond to : r = 1.25M (solid red), r = 1.5M 
(long dashed green), r = 2M (short dashed blue). The left panel shows the initial, trajectory dependent phase, and the right 
panel shows the later phase. Each trajectory is tracked until the shell crosses it. 



Using Eqn. (|28p the flux can be expressed as a function of time T, for an observer with r = const. The structure 
of the function S shows that the flux in the interior will remain zero as long as the shell trajectory satisfies Vs oc Vs, 
which is the case until t — 0. As a result, observers will see no flux as long as they are to the left of the ingoing ray 
which enters at r = 0. (See the Penrose diagram Fig. [31 This is in fact true outside the shell as well.) This is also 
confirmed by Fig. [6] which shows the flux as a function of time for three r = const observers. We see that after an 
initial complicated phase which depends on the trajectory details, the flux becomes positive and steadily increases in 
magnitude until the time the shell crosses the observer. This situation holds regardless of whether the shell crosses 
the observer before or after forming the event horizon, although the flux at shell crossing is larger for observers who 
are closer to the origin. 



B. Asymptotic trajectory 

In this section we discuss a class of trajectories first described in Ref. [Tsj . which are timelike throughout and 
asymptotically approach a final radius larger than the Schwarzschild radius of the shell. Classically the interior 
geometry is Minkowski spacetime and the exterior is Schwarzschild spacetime, with no event horizon. Although such 
a trajectory was discussed in detail in Ref. p^ . there are two reasons for returning to it here. Firstly, as mentioned 
earlier, the treatment in Ref. [l3l | was mainly analytical and hence the inversions required to calculate G{u) could 
not be performed. This can now be done in our current numerical setup, and we will display results that corroborate 
those in Ref. jl^ . Secondly, a numerical analysis allows us to further explore the semiclassical environment of this 
trajectory. Specifically, we will show below that such trajectories will generically give rise to a pulse of radiation which 
travels inwards, reflects at the origin and eventually reaches an external observer at late times. Mathematically, this 
new feature arises solely from the functional inversions at the shell boundary, and was (not surprisingly) missed in 
Ref. dl. 

The asymptotic trajectory is parametrized by Eddington retarded time u (which is a good parameter everywhere 
since there is no horizon), and is described by the following prescription, borrowed from Ref. [1^] : 

Ui{u)^e+(^(^l-^y -e^ e-"("' ee e + Ae~^^-^ , (31) 

where, as before, the prime is a derivative with respect to the argument and we are working with the non- 
dimensionalised form of the coordinates. Here < e < (1 — I/RqY^^ is a fixed constant, 

a{u) = ^ / duh{u) , (32) 
2 Jo 
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FIG. 7: The normalized exterior flux as a function of u, for the timelike asymptotic trajectory with Tio = 3M, e = 10"* and 
2MA = 0.25. The left panel shows the initial phase which shows a thermal behaviour followed by an exponential decay. The 
right panel shows the flux for the same conflguration at a much later time. The second, "pseudo-thermal" phase can be traced 
back to a pulse generated at an inflection point of the trajectory, as discussed in the text. 



where 0{u) is the Heaviside step function and h{u) is chosen to have the following asymptotic behaviour 



h{u) 
h{u 



0,M < 

■ oo) ^ 1 



h'{u) = 0,u < 0, 
h' {u — > oo) — > , 



(33) 



and we require the asymptotic values for h and h' to be achieved exponentially fast, with a time scale determined by 
M. For this paper, we choose h{u) = 6'(u) tanh^(Au) for some constant A. This is slightly different from the choice 
in Ref. [l^; this change does not alter the qualitative behaviour of the trajectory and was made to improve the 
numerical behaviour of various subsequent functions. The trajectory radius Rs satisfies the equation 



2i?l(l -U',) = Uf-[l- — 



1 



(34) 



so that Rs starts at a constant radius Rq, falls inwards and asymptotically approaches a final radius Rf = 1/(1 — e^). 
The Penrose diagram for such a trajectory was shown in Ref. (see their Fig. 3). 

We again consider two classes of observers : exterior observers at some constant radius r > i?o, and interior observers 
at some constant radius r < Rf. For brevity we ignore observers who are initially inside the shell and eventually 
come into the Schwarzschild region, since such observers will not see any effects that are absent for the previous two 
classes. For the exterior observers, the normalized fiux is easily seen to be identical in expression to the one calculated 
for the previous trajectory 



dG 



du 



where 



(35) 



(36) 



which is familiar from earlier, with r being replaced by u. The expression for Fu[dG/du] is somewhat simpler in terms 
of u. 



F„, 



dG 



du 



(37) 



Although these expressions may appear different from Eqn. (27) of Ref. [lj|, the two sets are in fact equivalent. The 
specific splitting used in Ref. [131 was conceptually easier to tackle, while the expressions above are easier to handle 
numerically. 

Since the dependence of the flux on the model parameters was discussed in detail in Ref. , here we will simply 
display results for a choice of parameter values that enhances the effects we want to focus on. The effects themselves 



12 



Flux 



1.x 10° 



5.x lO' - 



-5.x lO' - 



5.0 



(a) 



5.5 



6.0 



T/2M 



Flux 



1.x 10* 
8.x lO' 


: 




■ 






6.x lO' 


; 










4.x lO' 












2.x lO' 








I , , , 




2.x lO' 


5.14 


1 


5.16 


5.17 



T/2M 



(b) 



FIG. 8: The normalized interior flux as a function of Minkowski time T, for the timelike asymptotic trajectory with Ro = 3M, 
e = 10"* and 2MA = 0.25. The left panel shows the flux seen by three observers at constant r : r — OAM (solid red), r = 0.8M 
(long dashed black) and r — 1.6M (short dashed blue). The right panel shows details of the pulse (shown on the outgoing leg 
for the observer with r = OAM). 



are generic and do not depend on the specific parameter choices we make. Figure [7] shows the exterior flux as a 
function of u, for a trajectory with Rq = 3M, e = 10"'* and 2MX = 0.25. As expected for such a trajectory, we see 
an initial phase in the flux which closely resembles that in the conventional trajectory : the asymptotic trajectory 
initially behaves as if it were heading towards horizon formation. Eventually however, the effects of the nonzero e 
are felt and the flux exponentially decays to zero. This part is captured in the left panel of the figure, and was also 
described in Ref. [l3|- The right panel shows the flux for the same trajectory, but at much later times. We see 
a second pseudo-thermal phase (more on this nomenclature later), which rises from zero, reaches and stays at the 
Hawking value, and then also exponentially decays. The origin of this second phase will become clear soon, when we 
discuss the flux seen by the interior observers. 

For the interior observers, the flux calculation again proceeds identically to that for the horizon-forming trajectory. 
In fact, it is easily verified that the function S defined in Eqn. (j30p is parametrization independent, so that the 
normalized interior (outgoing) flux for the present case is given by 

^ = S{ui) - S{u2) , (38) 
1Q27r 

= - Pu[V:{u)]) , (39) 

with ui and U2 defined in terms of (T, r) via 

Vsiui)=T-r ; Vs{u2) = T + r . (40) 

As before, the flux will be nonzero only on the right of the ingoing ray that enters at u = 0, since in the region to the 
left of this ray we have (x Vg- Figure [5] shows the flux seen by three observers at three different constant radii r, for 
the same configuration as in Fig. [71 Each observer sees a strong incoming pulse of radiation and an outgoing pulse at 
a later time. By studying the time intervals between the pulses seen by various observers, it is possible to show that 
there is only one pulse which starts at the shell at a certain time, travels inwards and gets reflected at the origin, and 
is then seen by the observers again on its way out. This is also consistent with the fact that the second pulse seen by 
any interior observer is the precise mirror image of the flrst (although this is not obvious from the flgure). The right 
panel in Fig. [H] shows the detailed structure of the pulse. 

This discussion also suggests that this pulse should be detected by exterior observers as well, and in fact one flnds 
that the second pseudo-thermal phase referred to earlier, is precisely the remnant of this pulse after it exits the 
shell. Further, the origin of this pulse can be traced to the time interval in which the behaviour of Us changes from 
~ (1 — l/RoY^'^u to Us eu (see Eqn. ([?T|) ). To visually demonstrate this we deflne the time Wchangc such that 
the slope C/^ (^change) — ee, i.e. one e-folding away from its final asymptotic value. In Fig. [51 the vertical dashed line 
keeps track of the ray which enters at Uchangc- This figure uses the same trajectory configuration as in Figs. [3 and 
m except for the value of e, which is now e = 10""^. Figure [9l when compared with Figs. [71 and [51 also suggests that 
the pulse height in the interior scales approximately as e~^, which also follows from considering the structure of S, 
and that the time at which the second phase occurs in the exterior, scales roughly as e~^. This is also supported by 
a more detailed numerical analysis (not shown). 
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FIG. 9: The behaviour of the pulse which arises in the interval when Us changes its behaviour (panel (a)), for the asymptotic 
trajectory configuration Ro — 3M, e = 10"'^ and 2MA = 0.25. The vertical dashed line in panel (a) marks u = Uchango, and the 
vertical lines in the other panels track the ray TZ that enters the shell at « = Uchangc- Panel (b) : the interior normalized flux 
for the observer at r = OAM. Vertical lines mark the times when this observer meets the ray TZ, before and after reflection. 
Panel (c) : The initial phase of the exterior flux. Vertical line marks the ray TZ on the way in. Panel (d) : The late phase of 
the exterior flux. Vertical line marks the ray TZ on the way out. 



One significant aspect of the pulse of radiation is that it is very sharply peaked in the interior of the shell, and 
becomes diffuse once it exits to the exterior. In fact one can check that the height of the pulse in the exterior never 
exceeds the Hawking value, instead the width increases with decreasing e, and the flux in the exterior pulse spends 
a longer time at the Hawking value. This is why we referred to this phase as "pseudo-thermal" earlier. It will be 
very interesting to see the signatures left by this pulse in the spectrum of particles in various regions. This however 
requires a calculation of the Bogolubov coefficients, which is complicated by the presence of numerical inversions, and 
is beyond the scope of this work. Another tantalizing prospect is that such a pulse could possibly lead to astrophysical 
signatures in objects which undergo phases of halted collapse without forming a horizon. This would, however, require 
some kind of amplification mechanism, since the total amount of gravitational energy released in the pulse appears 
to be very small (while the pulse height can become large, the pulse width decreases with decreasing e). 



IV. DISCUSSION 



We have explored in detail the scmiclassical environment of a collapsing shell of matter. Our main objective was to 
study the scmiclassical flux observed by a variety of observers, both inside and outside the shell. We have studied 
two types of trajectories : those which form a horizon and those which do not. 

For the horizon forming trajectories, we reproduced the well known result of an asymptotically constant flux seen 
by observers outside the horizon. Additionally we demonstrated that freely falling observers measure an outgoing 
flux which generically increases without bound as they approach the singularity. Our calculation assumed zero 
backreaction, and it is possible that the infinity we see in our flux can somehow be cured if backreaction were correctly 
accounted for throughout. It is hard to guess how exactly this might happen, since a self-consistent calculation of a 
black hole geometry in the presence of backreaction at late stages of evaporation, does not exist. It is possible that 
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the resolution of this problem is linked to the more general problem of the singularity itself. 

We also saw, at least in the timelike case, that inertial observers inside the shell also see a nonzero outgoing flux, 
which steadily increases until the time they exit the shell. In the context of backreaction, this has consequences for 
all arguments built on the assumption that the interior geometry of the shell remains that of Minkowski spacetime 
throughout. This assumption was made for example in Ref. [isj while arguing that the semiclassical backreaction can 
at best delay the formation of the event horizon by a short amount of time. While it is unlikely that including a flux 
of the same order of magnitude as the Hawking flux would overturn the conclusions of the arguments in Ref. fisj , it 
is nevertheless important to rigorously account for this additional flux and build a complete picture. Speciflcally, it 
is interesting to ask what the self-consistent geometry in the interior will be, after accounting for backreaction. We 
will return to this question in future work. 

In the case of the asymptotic trajectory which halts collapse before horizon formation, we corroborated the analysis 
of Ref. by a full numerical calculation, and extended the results to include the flux seen by interior observers. In 
the process we uncovered an interesting, generic feature of such trajectories, in the form of a pulse of radiation which 
is sharply peaked in the interior but appears to thermalize as it emerges from the shell. It remains to be seen whether 
this feature might lead to astrophysically interesting signatures. 

For both types of trajectories discussed here, it will be interesting to determine the nature of the particle spectrum, 
for which we need to calculate the Bogolubov coefficients in the respective geometries. This is work in progress. We 
expect our current results to be useful in addressing the broader issues mentioned in the Introduction, specifically 
questions related to the interpretation of thermodynamic variables such as entropy, for a variety of classes of observers, 
and in studying the role and nature of the semiclassical backreaction. 
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